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Abstract 

When  a  wave  of  wavelength  X  is  incident  upon  an  opaque  object  of  typical 
dimension  a,  a  shadow  is  formed  in  the  geometric  optics  limit  X/a  =0.   If  X/a 
is  small  and  not  zero,  the  shadow  boundary  is  displaced  slightly  from  the  geo- 
metrical shadow  boundary.   It  is  known  from  the  work  of  Artmann  that  the  amount 

of  the  shift  in  the  case  of  circular  cylinders  is  asymptotically  equal  to 

2  1  / "*> 
a(X  a)    .   The  proportionality  constant  a  may  be  positive  or  negative,  de- 
pending on  the  boundary  conditions  satisfied  by  the  field.  A  similar  result 
was  obtained  for  parabolic  cylinders  by  Rice. 

Because  of  the  discrepancy  in  the  proportionality  constant  for  hard  cylinders 
that  exists  between  the  works  of  Rice  and  Artmann,  already  pointed  out  by  Wait 
and  Conda,  we  have  redetermined  the  shift  for  the  circular  cylinder  and  found 
that  a  was  the  same  as  for  the  parabolic  cylinder  as  determined  by  Rice.   In 
addition,  the  shift  has  been  determined  for  circular  cylinders  for  fields  which 
satisfy  an  impedance  boundary  condition. 

The  agreement  of  the  results  for  circular  cylinders  and  parabolic  cylinders 
can  be  understood  in  terms  of  the  local  nature  of  shadow  boundary  rays.   On  this 

basis  we  conjecture  that  the  shift  of  the  shadow  boundary  cast  by  any  two  or  three 

2l/5 
dimensional  object  when  X/a  is  small  is  asymptotically  equal  to  a(X  a)    along 

each  boundary  ray.  The  constant  a  depends  upon  the  impedance  of  the  object  at  the 
point  of  tangency  and  is  given  by  our  result  for  the  circular  cylinder.  This  con- 
jecture is  consistent  for  waves  at  oblique  incidence. 

We  have  also  determined  the  scattering  cross  section  per  unit  length  of  a 
circular  cylinder  for  X/a  small  and  for  an  impedance  boundary  condition.   We  find 
that  the  deviation  from  the  geometrical  optics  cross  section  is  proportional  to 
the  shift  of  the  shadow  boundary.  This,  together  with  the  local  character  of  the 


of  the  shift,  leads  us  to  conjecture  that  the  total  scattering  cross  section  for 
any  two  or  three  dimensional  object  when  X/a  is  small  is  asymptotically  equal  to 
the  integral  around  a  cross  section  of  the  shadow  boundary  of  a  multiple  of  the 
shift.   This  conjecture  is  checked  by  applying  it  to  spheres  and  showing  that  it 
yields  the  known  exact  results.   Finally,  we  show  how  similar  results  and  con- 
jectures may  be  applied  to  the  case  of  incident  electromagnetic  waves. 
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1.   Introduction. 

When  a  wave  Is  incident  upon  an  opaque  object  of  typical  dimension  a, 
vhich  is  large  compared  to  the  incident  wavelength  X,  a  shadow  is  formed. 
In  the  geometrical  optics  limit,  X/a  =  0,  the  shadow  is  bounded  by  those 
optical  rays  which  are  tangent  to  the  object.  However  when  X/a  is  small 
and  not  zero,  the  shadow  boundary  is  parallel  to,  but  displaced  slightly 
from,  the  geometrical  shadow  boundary.  ThiB  fact  was  first  discovered 
theoretically  by  K.  Artmann  in  1950^-1.   He  considered  a  scalar  plane  wave 

normally  incident  upon  a  circular  cylinder  of  radius  a  and  showed  that  the 

2  1/5 
shift  was  asymptotically  equal  to  <x(X  a)  '    .  The  proportionality  constant 

a  was  found  to  be  positive  for  a  soft  cylinder,  on  which  the  field  vanishes, 

and  negative  for  a  hard  one,  on  which  its  normal  derivative  vanishes.  Thus 

in  the  former  case  the  shadow  was  enlarged  while  in  the  latter  case  it  was 

r2"i 

diminished.  A  similar  result  was  later  obtained  by  S.O.  Rice  L  J  for  the 
shadow  cast  by  a  parabolic  cylinder.   In  this  case  a  was  the  radius  of 
curvature  of  the  cylinder  at  the  point  of  tangency  of  the  geometrical  shadow 
boundary. 

Since  fields  propagate  along  rays  when  X/a  is  small,  we  then  expect  each 
ray  forming  a  shadow  boundary  to  be  shifted  indepently  of  other  such  rays . 
In  other  words,  we  expect  the  shadow  boundary  shift  to  be  a  local  effect.   If 
so,  the  shift  of  each  shadow  boundary  ray  should  depend  only  upon  the  curvature 
of  the  object  at  the  point  where  the  ray  is  tangent  to  it.   This  expectation 
is  partly  confirmed  by  the  fact  that  Artmann  obtained  nearly  the  same  value  for 
the  constant  a  for  a  soft  circular  cylinder  as  Rice  did  for  a  soft  parabolic 
cylinder.   However  they  obtained  unequal  values  for  hard  cylinders.   Since  this 
result  is  contrary  to  our  expectation,  and  since  Rice  checked  his  work  without 


finding  an  error,  we  redetermined  the  shift  for  the  circular  cylinder  and 
found  that  a  was  the  same  as  for  the  parabolic  cylinder  in  both  cases.  Thus 
the  expectation  that  the  shift  is  a  local  effect  is  confirmed.   In  addition 
we  have  determined  the  shift  for  circular  cylinders  on  which  the  field  satisfies 
an  impedance  boundary  condition. 

This  confirmation  leads  us  to  conjecture  that  the  shift  of  the  boundary 

of  the  shadow  cast  by  any  two  or  three  dimensional  object  when  X/a  is  small 

2  '13 

is  asymptotically  equal  to  a(X  a)  along  each  boundary  ray  (see  equation  (27)). 

Here  a  denotes  the  radius  of  curvature  of  the  surface  of  the  object  at  the 
point  of  tangency  of  the  ray,  in  a  plane  containing  the  ray  and  the  normal  to 
the  surface.  The  constant  a  depends  upon  the  impedance  of  the  object  at  the 
point  of  tangency  and  is  given  by  our  result  for  the  circular  cylinder.  To 
check  our  conjecture  we  show  that  it  is  consistent  for  a  plane  wave  obliquely 
incident  upon  a  cylinder  of  arbitrary  cross  section.  This  means  that  it  is  true 
for  oblique  incidence  if  it  Is  true  for  normal  incidence.   Since  it  is  true  for 
normal  incidence  in  the  cases  of  a  circular  cylinder  of  any  impedance  and  of 
soft  and  hard  parabolic  cylinders,  it  is  also  true  for  oblique  incidence  in 
those  cases. 

Next  we  determine  for  X/a  small,  the  total  scattering  cross  section  per 
unit  length  of  a  circular  cylinder  on  which  the  field  satisfies  an  impedance 
boundary  condition.  We  find  that  the  deviation  from  the  geometrical  optics 
cross  section  is  proportional  to  the  shift  of  the  shadow  boundary.   This 
proportionality,  and  the  local  character  of  the  shift,  suggest  a  conjecture 
for  the  total  scattering  cross  section  of  any  two  or  three  dimensional  object 
when  X/a  is  small.   It  is  that  the  deviation  from  the  geometrical  optics  cross 
section  Is  asymptotically  equal  to  the  integral,  around  a  cross  section  of  the 
shadow  boundary,  of  a  multiple  of  the  shift.  The  multiplier  is  just  the  proportionality 
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factor  found  for  the  circular  cylinder.   This  conjecture  is  checked  by 
applying  it  to  hard  and  soft  spheres  and  showing  that  it  yields  the  known 
exact  results. 

Finally  we  consider  electromagnetic  waves  incident  upon  perfectly  con- 
ducting objects.   We  determine  the  shadow  boundary  shift  and  the  scattering 
cross  section  for  a  circular  cylinder.  We  obtain  results  similar  to  those 
we  obtained  in  the  scalar  cases  and  make  similar  conjectures  for  arbitrary 
objects. 

The  field  near  the  boundary  of  the  shadow  of  a  circular  cylinder  has 
also  been  determined  by  J.R.  Wait  and  A.M.  Conda  [3]  when  the  field  satisfies 
an  impedance  boundary  condition  on  the  cylinder.   In  the  cases  Z=0  and  Z=co 
they  have  shown  that  the  value  of  the  field  on  the  geometrical  optics  boundary 
agrees  with  that  found  by  S.O.  Rice  [2]  on  the  geometrical  optics  boundary  of 
the  shadow  of  a  parabolic  cylinder.  They  have  also  pointed  out  that  the 
corresponding  results  of   K.  Artmann  [l]  differ  from  theirs  and  Rice's,  and 
attribute  this  difference  to  an  inaccuracv  in  his  calculation. 
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Figure  1   -  A  plane  wave  e    is  incident  from  the  left  upon  a  circular 
cylinder  of  radius  a. 


Figure  2   -  The  shadow  cast  by  a  circular  cylinder  illuminated  from  the  left 
by  a  plane  wave.  The  dashed  lines  are  the  geometrical  shadow 
boundaries.   The  solid  lines  bound  the  transition  zones  between 

the  illuminated  and  shaded  regions .   The  width  of  each  zone  is 

l/2 
approximately  2(2x/k)    at  a  distance  x  behind  the  cylinder. 

The  two  zones  merge  at  the  (approximate)  distance  ka  /2  which 

determines  the  length  of  the  shadow. 


2.  The  scattered  field. 


We  wish  to  determine  the  field  u(r,9)  which  results  when  a  plane  wave 

ikx  / 

e    is  incident  from  the  left  upon  a  circular  cylinder  of  radius  a.   (See 

Figure  1.)  The  field  satisfies  the  reduced  wave  equation 


(A  +  k  )u  =  0,    r  5  a 


(1) 


On  the  cylinder  it  satisfies  the  impedance  boundary  condition 


.am 


fiu  =  (-st-  +  iZu)  =  0  at  r  =  a 


(2) 


Here  Z  is  a  constant  called  the  impedance  of  the  cylinder  surface  and  ft  is  the 

ikx 
boundary  operator.   The  scattered  field  u  (r,9),  defined  by  u  =  e    +  u  , 


must  satisfy  the  radiation  condition 
ou 


lim 


vdr      s' 


(5) 


Equations  (l)  -  (3)  have  a  unique  solution  u  for  which  u  is  given  by 


«B(r,e)  "  -  |  5~   e  Hn  V) 


nH^2)(ka)" 


HH^^ka) 


(*) 


We  can  rewrite  (4)  in  another  form  by  introducing  F(r,9),  the  Fourier 
transform  of  the  summand,  defined  by 


F(r 


..)--*   / 


iv(f  +  9)  ,.» 

v   v   ' 


1  + 


^2)(ka)- 


V1}(ka)_ 


(5) 


In  terms  of  F,  the  Poisson  summation  formula  yields  for  (k)} 


U  (r,9)  =  YZ       F(r,2:tm+9) 
m=  -oo 


(6) 
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We  now  rewrite  the  terms  in  (6)  for  which  m  is  negative  hy  replacing  m 

by  -m'.  Then  upon  changing  v  to  -v '  in  (5)  and  making  use  of  the  relations 

between  Hankel  functions  of  positive  and  negative  orders,  we  finally  obtain 


u  (r,0)  =  G(r,9)  +  G(r,-6) 


m=l 


[F(r,2rfm+0)  +   F(r,2jtm-e)] 


Here  G(r,9)  is  defined  by 


G(r,9)  =  - 


H^;(kr) 


V2)(ka)l 
h(v)  +  ^rj dv 


(7) 


(8) 


In  (8)  h(v)  is  the  unit  step  function  defined  by 
h(v)  =  0,  v  <  0;  h(v)  =1,  v  >  0 


(9) 


From  u  we  can  compute  <f,   the  scattering  cross  -section  per  unit  length 
of  cylinder.   If  Z  is  real  the  cross  section  theorem  yields 


=  lim   x.p.feV^-i^^A)   (r,0)] 
r  ->  00     L  k  s    J 


R.P., 


fiH(2)(ka)~ 
ti(v)  +-7^ 


dv  +±B.?.TZ 

m=l 


1  + 


n^2)(kay 


SP(ka) 


dv 


(10) 


When  Z  is  complex  C  denotes  the  extinction  cross  section  which  is  the  sum  of 

the  scattering  and  absorption  cross  sections. 

3.   Asymptotic  form  of  u  near  the  shadow  boundary. 

The  term  F(2toh+€))   in  (7)  represents  a  diffracted  or  creeping  wave 
which  has  traveled  m  times  around  the  cylinder  in  the  direction  of  increasing 
9.   Similarly  F(2joa-9)    represents  a  wave  which  has  traveled  m  times  around 


in  the  direction  of  decreasing  9.  These  diffracted  waves  decay  exponentially 
with  increasing  angle  around  the  cylinder  with  a  decay  rate  proportional  to 
(ka)  '  .  Therefore  for  large  values  of  ka,  which  means  for  a  large  compared 
to  the  wavelength  X  =  2*/k,  all  the  terms  in  the  sum  in  (7)  are  exponentially 
small  compared  to  G.   Consequently  in  evaluating  ug  for  ka  large,  we  need  only 
consider  the  terms  G(r,9)  and  G(r,-9). 


We  shall  find  it  convenient  to  split  the  integral  in  (8)  into  two 
terms  G1(r,9)  and  G  (r,9)  defined  by 


";;  iv(f  +  9)  m 

^,9)  «  -|  J     e        H^Vrjdv  , 


(11) 


G0(r,6)  = 


iv(|  +  9)  .    ,     n^2)(ka) 


OH^fka) 
v   v   ' 


1  T     ±V(2   +  0)  (1) 

I  J     •   2     ^(kr) 


nH(2)(ka)" 
v  v   ' 


nH(l)(ka) 
v  v  '- 


(12) 


Then 


G(r,9)  =  Gx(r,9)  +  G2(r,e) 


(15) 


To  evaluate  G,  asymptotically  for  large  ka  and  r  >  a  we  may  replace  IT  (kr) 
by  its  Debye  asymptotic  form  which  is  valid  when  kr  is  large  and  kr  >  v .  We 
also  introduce  the  integration  variable  t  =  sin_1(v/kr).  Then  (ll)  becomes 


sin"    (a/r) 
G;L(r,e)~  -e"ln/1+(g)l/2      f  (cos   T)l/2exp|itr[cos   t  +  (e+T)sin  x]|  dx 

o 

(1*0 

Let  us  now  evaluate  G(r,-9)  on  and  near  the  upper  boundary  of  the  shadow 
cast  by  the  cylinder.  Since  ka  is  large  and  r  >  a,  kr  is  also  large  so  we  may 
employ  the  method  of  stationary  phase  to  evaluate  G  (r,-e).   The  stationary  phase 
point  is  at  t=9  or,  in  terms  of  the  cartesian  coordinate  y,  at  T=sin  (y/r) . 
Near  the  shadow  boundary  y=a  this  stationary  point  is  near  the  upper  limit  of 
integration.  This  must  be  taken  into  account  in  evaluating  the  integral 
asymptotically.   By  expanding  the  phase  up  to  quadratic  terms  around  the 
stationary  point  we  obtain 

sin"  (a/r) 
Gi(r,.6)  ^..-iAg)l/2-(f)l/2  ."*  f  elk^-6>  '2   d.        (15) 

O 

l/2 

Now  we  introduce  t  =  (f-e)(kx/2)  /   as  integration  variable  and  (lj)   becomes 


_  ,     M        i(kx-«A)  -1/2    r 

G.(r,-e)  -v  -e  N    '  '«  '    / 


(kx^J^Qsin-^a/rj-e] 


(16) 


-(kx^-^e 


If  the  lower  limit  of  the  integral  in  (l6)  is  large  in  absolute  value  and 
the  upper  limit  small  in  absolute  value,  (l6)  yields 

~  -  V  |I  ♦  e-^Cg^a^)]  (17) 

The  conditions  under  which  (17)  holds,  which  were  imposed  upon  the  integration 
limits  in  (l6),  are 


(|£)l/2  ten-l(a/=0  »  1  ,  (18) 

(|;)1/2  |a-y|  «  1     .  (19) 

Equation  (l8)  is  equivalent  to  x  «  ka  /2,  which  means  that  x  must  be  small 

compared  to  the  length  of  the  shadow  ka  /2.   Equation  (l?)  requires  that 

the  distance  |a-y|  of  the  point  from  the  geometrical  shadow  boundary  y=a 

be  small  compared  to  the  half  width  (2x/k)  '   of  the  shadow  boundary  zone.   (See 

Figure  2). 

To  evaluate  G?(r,-©)  asymptotically  for  large  ka  we  first  insert  into 
the  integral  the Debye  asymptotic  form  for  the  Hankel  functions.  This  asymptotic 
form  is  valid  only  when  the  order  and  argument  are  unequal.   Therefore  it  is  not 
valid  for  one  of  the  Hankel  functions  when  v=kr  and  for  the  other  two  when 
v=ka.   However  the  neighborhoods  of  these  two  points  contribute  asymptotically 
negligible  amounts  to  the  integral  if  the  point  r,0  does  not  lie  on  or  near  the 
upper  boundary  of  the  shadow.  This  can  be  seen  by  noting  that  at  all  other  points 
in  the  illuminated  region  for  which  0  <  0  <  it,  the  integral  can  be  evaluated 
asymptotically  by  the  method  of  stationary  phase.   The  stationary  point,  whose 
neighborhood  yields  the  entire  asymptotic  contribution  to  the  integral,  is  away 
from  v=kr  and  v=ka.   In  the  lower  half  of  the  (r,0)  plane,  the  entire  Integral 
is  asymptotically  negligible  compared  to  G„(r,0).   However  when  the  point  r, 0 
approaches  the  upper  boundary  of  the  shadow,   the  stationary  point  approaches 
v=ka.  Therefore  we  conclude  that  for  points  on  and  near  this  boundary,  the 
entire  asymptotic  contribution  to  the  integral  comes  from  the  neighborhood  of 
v=ka.   But  then  the  Debye  asymptotic  forms  of  H^   (ka)  and  IT   (ka)  cannot  be 
used  since  they  are  not  valid  at  v=ka.   Instead  we  must  use  their  uniform  asymptotic 
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representations  in  terms  of  the  Airy  function  A.  For  IT   (kr)  we  may  still 
use  the  simpler  Debye  asymptotic  form  provided  that  r  >  a,  and  expand  it 
around  v=ka.   In  the  special  case  of  forward  scattering,  for  which  0=0  and 
r  »  a,  it  has  been  proved  that  this  procedure  does  yield  the  asymptotic  form 
of  the  integral  IAI  . 

In  order  to  carry  out  the  procedure  just  described  for  points  on  and  near 
the  upper  shadow  boundary  we  first  consider  the  Debye  asymptotic  form  of 
IP   (kr) .   In  it  we  replace  v  by  the  new  variable  q  by  means  of  the  equation 


v  =  ka  +  (^)l/5q  (20) 

iv(|  -  e) 

Then  we  multiply  it  by  e         and  simplify  the  product  by  considering 
y-a  to  be  small  and  x  to  be  positive.   In  this  way  we  obtain 

Now  we  insert  (21)  and  the  uniform  asymptotic  forms  of  the  other  Hankel 
functions  into  (12)  and  introduce  q  as  the  integration  variable.  Then  we  obtain 
in  the  neighborhood  of  the  upper  shadow  boundary 

G2(r,-9)  ~-(2«kxrl/2e   "  k   (ka^O^ka)1/^  .  (22) 

Here  C  is  defined  by 


^h]  -  e^fibw&g* ♦  6-1/5/{i  +  f [q,(M1/3z] 


OD 

dq  (23) 


where 


f  [q,(ka)l/3SQ  -  e^  A'(A)  -  ie"l5(ka/6)1^ZA(e1^)   _  (A) 

-I-      i-  -i— 

A*(e   3q)  -  ie  5(ka/6)1/5ZA(e   5q) 


The  Airy  function  in  (2k)    is  defined  by 
oo 


A(q)  =  |  J        exp[-i(|5-qi)]d|  (25) 


For  points  at  which  the  conditions  (l8)  and  (19)  are  satisfied,  G2(r,0) 
is  asymptotically  negligible  compared  with  G?(r,-0).  Physically  this  is  so 
because  Gp(r, 0)  describes  the  reflected  wave  in  the  region  0  >  9  >  -«  and 
describes  a  diffracted  wave  in  the  region  0  <  0  <  it.  The  diffracted  wave 
decays  exponentially  with  increasing  0  and  it  is  exponentially  small  compared 
to  Gp(r,-0)  near  the  upper  shadow  boundary,  provided  (l8)  and  (19)  are  satisfied. 
Thus  when  these  conditions  are  satisfied  the  field  on  and  near  the  upper  shadow 
boundary,  u  =  e    +  u  ,  is  from  (7),  (13),  (17)  and  (22) 

u„  e^  .  (2rtkx)-1/2ei^kX-Jt/4)  {k(a-y)+(ka)1/5c[(ka)1/5z]}  .         (26) 

k.      The  displacement  of  the  shadow  boundary. 

From  (26)  we  see  that  at  the  geometrical  shadow  boundary  y=a,  u  approaches 
e   /2  as  ka  becomes  infinite.  Thus  the  field  on  the  geometrical  shadow  boundary 
approaches  one  half  the  incident  field  in  the  geometrical  optics  limit.  This 
suggests  that  we  define  the  shadow  boundary  for  finite  values  of  ka  as  the 
curve  on  which  |u|  =  l/2.   By  applying  this  definition  to  (26)  we  find  that 
the  shadow  boundary  is  then  the  straight  line  y  =  a+s  where  3  is  given  by 

s  =  (CR+Cl)(ka)-2/5a  (27) 

Here  CR  and  C  are  the  real  and  imaginary  parts  of  C[(ka)  '  z]  .  Thus  the 
shadow  boundary  is  displaced  or  shifted  from  the  geometrical  boundary  y=a  by 
the  amount  s  which  we  shall  call  the  shift. 


-  11  - 


In  1950  K.  Artmann'--i  showed  that  for  Z=0  and  Z=ao  the  boundary  of  the 
shadow  of  a  circular  cylinder  was  shifted  by  the  amounts 

s  =  -.20(2rtka)-2/5a  -  -  .68(ka)_2/5a      z  =  O^^1   „  0] 


(28) 


s  =  .39(2jrka)"2Ma  =  1.33(ka) ~2' 5a       z  =  oo'[u(a,9)  =  6]  (29) 

In  1954  S.O.  Rice L J  showed  that  for  Z=0  and  Z=co  the  boundary  of  the  shadow 
cast  by  a  parabolic  cylinder  was  also  shifted.   In  terms  of  the  radius  of 
curvature,  a,  of  the  parabola  at  the  point  of  tangency  of  the  geometrical 
shadow,  the  shift  was  found  to  be 

s  =  -.346(2*ka)"2/5a  =  -l.ie(ka)"2/3a     Z  =  0  (50) 

s  =  .399(2rtka)"2/3a=  1.56(ka)~2/3a     Z  =  00  (31) 

The  process  of  shadow  boundary  formation  is  local  in  character.  Therefore 

it  is  to  be  expected  that  the  leading  term  in  the  expression  for  the  shift 

should  depend  only  upon  the  radius  of  curvature  of  the  shadow  forming  cylinder 

at  the  point  of  tangency  of  the  geometrical  shadow  boundary.   The  results  (29) 

and  (31)  tend  to  confirm  this  expectation  while  (28)  and  (30)  seem  to  contradict 

it.   Let  us  compare  these  results  with  our  result  (27).  When  Z=0  and  Z=co , 

S.I.  Rubinow  and  T.T.  Wu^-J  and  T.T.  Wu  ^   have  shown  that  C  =  >/3  CR  and 

have  computed  (^(0)  =  -0.432120  and  (1,(00)  =  0. 498077-  Using  these  results 
R  " 

in  (27)  we  obtain  for  the  shifts  in  these  cases 

s  =  -1.180574   (kaf  2/3a  Z  =  o["M^9)   =  0]  (32) 

s  =  1.360771  (ka)"2/5a  Z  =  co  [u(a,e)   =  0]  (35) 


!:■ 


We  see  that  our  results  (32)  and  (33)  for  the  circular  cylinder  agree  with 
those  for  the  parabolic  cylinder  given  in  (30)  and  (3l).  Thus  the  discrepancy 
was  due  to  an  inaccurate  determination  by  Artmann  of  the  constant  in  (28). 

This  agreement  confirms  our  expectation.   It  also  leads  us  to  conjecture 
that  the  shift  of  the  boundary  of  the  shadow  cast  by  any  object  is  given  by 
(27)/  as  described  in  the  Introduction.  To  check  it,  let  us  consider  a  plane 
wave  incident  obliquely  upon  a  cylinder  of  arbitrary  cross  section.   Let 
rt/2  -  7   denote  the  angle  between  the  incident  rays  and  the  generators  of  the 
cylinder.   Let  a'  be  the  radius  of  curvature  of  the  normal  cross  section  of 
the  cylinder  at  one  of  its  two  points  of  tangency  with  the  shadow  boundary. 

Then  the  radius  of  curvature  a  of  the  cylinder  surface  in  the  plane  containing 

2 
the  incident  tangent  ray  and  the  normal  to  the  cylinder,  is  a  =  a' sec  7. 

Thus  (27)  yields  for  the  shift 

s  =  (C^Kka'r^VCsec  7f/3 

where 

C  =  c[(ka')l/5Z(sec  7)2/3l 

Another  expression  for  s  can  be  obtained  from  the  solution  v(x,y, z)  of  the 
boundary  value  problem  for  oblique  incidence.   It  is  known  that  v  can  be  related 
to  u(x,y),  the  solution  for  normal  incidence,  which  we  have  previously  considered. 
The  relation,  which  involves  replacing  k  by  k  cos  7   and  z  by  z  sec  7, Is 

v  =  erp[ikz  sin  7~]   u(x,y,k  cos  7,    z  sec  7). 

If  (27)  is  used  for  the  shift  of  the  shadow  boundary  of  u(x,y,k  cos  7 ,    2  sec  7), 
it  yields  exactly  the  result  given  above.  This  agreement  shows  that  our  con- 
jecture is  consistent,  i.e.,  it  Is  correct  for  oblique  incidence  if  It  Is  correct 
for  normal  Incidence.   In  those  cases  for  which  we  know  it  to  be  correct  at 
normal  incidence  (circular  and  parabolic  cylinders)  it  is  therefore  also  correct 
for  oblique  Incidence. 
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5.  The  scattering  cross  section  of  the  cylinder. 

When  ka  is  large,  the  cross  section  o"  is  asymptotically  given  by  the 
first  term  in  (10)  since  the  sum  is  relatively  negligible  for  the  reasons 
outlined  previously.  That  term  is  an  integral  which  can  be  evaluated  by 
splitting  it  into  Oi   and  cf?  just  as  the  integral  for  G(r,G)  was  split  in 
(11)  and  (12). 


QEK    '(ka) 
-oo    v  v  ' 


I 


OH^(ka) 

1  +   (Ti 

me,  (ka) 


(3M 


dv  >  (35) 


Then 


(36) 


We  note  that  ol  =  ka.   is  twice  the  diameter  of  the  cylinder,  so  it  is  just 
the  geometrical  optics  scattering  cross  section.   Thus  cf    is  the  correction 
to  it  due  to  the  finiteness  of  ka. 

The  integrals  in  (35)  can  be  evaluated  in  the  same  way  as  were  those  for 
G2(r,e)  in  (12).   The  result  is  similar  to  (22),  and  (35)  becomes 

<f2   =  £(ka)l/5CR  (37) 

Here  C  is  the  real  part  of  c[(ka)1'5z]  which  is  given  by  (23).   Now  (36)  becomes 

<T~  4a  f"l  +  CR(ka)"2/5l  (38) 

M 
For  Z=0  and  Z=oo  this  result  coincides  with  that  of  S.I.  Rubinow  and  T.T.  WuLJ. 

6.  Relation  of  the  shadow  boundary  shift  to  the  cross  section. 
Upon  comparing  (27)  with  (38)  we  find 


Ik 


<f~  k 


(39) 


Equation  (39)  shows  that  the  scattering  cross  section  is  the  same  as  the 

CR 
geometrical  optics  cross  section  of  a  cylinder  with  the  radius  a  +  - — —  s. 

CR+CI 

For  Z=0  and  Z=oo ,  CT.(CD+CT)"1  =  (l  +  V^)"1- 
K   n   1 

Let  us  now  consider  a  non-circular  cylinder  with  two  shadow  boundaries 
and  two  shifts  s  and  s?.  We  expect  each  shift  to  be  given  by  (27)  with  the 
appropriate  radius  of  curvature  a  or  a  of  the  cylinder  at  its  point  of 
tangency  with  the  corresponding  shadow  boundary.   Let  w  denote  the  width  of 
the  shadow  according  to  geometrical  optics.  Then  by  analogy  with  (39) >    and 
based  upon  the  physics  of  shadow  formation,  we  expect  the  scattering  cross 
section  of  the  cylinder  to  be 


(' 


:-2/'(CRa1/5)i  +  ,T2/'(C//5 


Q- 


(ho) 


In  case  Z=0  or  Z=co.  CL  is  independent  of  a,  and  (ko)    simplifies  slightly. 
This  simpler  formula  was  suggested  previously  on  the  basis  of  similar 
considerations  L  -•  . 

Next  we  consider  the  shadow  of  a  three  dimensional  object.  Along  each 
ray  of  the  geometrical  shadow  boundary  we  expect  a  shift  normal  to  this  boundary. 
The  shift  should  still  be  given  by  (27)  with  a  the  radius  of  curvature  of  the 
object's  surface  at  the  point  of  tangency  of  the  ray,  in  the  plane  containing  this 
ray  and  the  surface  normal.   Let  t  denote  arclength  along  a  normal  cross 
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section  of  the  shadow  boundary.  Then  if  A  is  the  cross  sectional  area  of  the 
geometrical  shadow  we  should  have 

C. 


5  rf)  C^dTl 


[A+^c^dT]  =2[A  +  k" 


This  expression  as  well  as  (ko)    should  also  apply  when  Z  is  not  constant. 

When  Z=0  or  Z=oo  and  C_  is  independent  of  a,  (kl)   becomes 
K 

0--2[A  +  CRk-2/3  <£   a^V].  (42) 


As  an  example  of  the  use  of  (k-l)   let  us  determine  <S  for  a  "body  of 
revolution  with  the  incident  plane  wave  traveling  along  its  axis.   Let  p 
denote  the  maximum  radius  of  the  body.  Then  (kl)   becomes 

C~  2*p2  +  WCDk-2/3al/5p  (k-3) 

H 

For  a  sphere, p=a, and  (kj)   reduces  to 

o"~  2jta2KL  +  2CR(ka)"2/5 

In  the  special  cases  Z=0  and  Z=oo  the  result  (kk)   for  the  sphere  agrees 
with  previously  known  results  [1»-,5j'  This  agreement  is  a  verification  of 
the  assumptions  on  which  (kl)    is  based. 


m 


1> 


7.   The  electromagnetic  case. 

Let  us  now  consider  an  elliptically  polarized  monochromatic  plane 
electromagnetic  wave  incident  normally  upon  a  perfectly  conducting  cylinder 
of  arbitrary 
is  given  by 

EinC  =  (k  cos  3  +  "J  e15  sin  p)eikx  (1+5) 

Here  "J  &nd  k  are  unit  vectors  along  the  y  and  z  axes,  respectively,  while  (3 

and  8  are  real  constants.   The  cylinder  generators  are  parallel  to  the  z  axis 

-icot  1 — 

and  the  time  factor  e     is  omitted.   The  propagation  constant  is  k  =  cu  sjei± 

where  e  and  \i   are  the  dielectric  constant  and  permeability  of  the  medium.  Then 

the  total  fields,  electric  and  magnetic, consisting  of  the  incident  and  scattered 

fields,  can  be  written  as 

E  =  ku   cos  3  -  ik~1ei6  sin  p  k  X  Vq  (k6) 


H  =  \J  -       (kuQ  e16  sin  p  +  Ik"1  cos  pkX  VuJ 

In  these  equations  u  and  u   are  the  fields  which  result  when  the  scalar 
^  o       CO 

ikx 
plane  wave  e    is  incident  upon  a  hard  or  soft  cylinder  with  the  same  cross 

section  as  the  given  one. 

On  and  near  the  upper  boundary  of  the  shadow  of  a  circular  cylinder,  u 

and  u   are  given  by  (26).   This  equation  shows  that  in  both  cases  Vu  ~  ikui, 

where  i  is  the  unit  vector  along  the  x  axis.  Upon  using  this  result  in  (46) 

and  (J+7)  we  obtain  for  the  fields  on  and  near  the  shadow  boundary 


(4?) 
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E  -  ku     cos  p  +  ju  e       sin  p  (48) 


H-"    \l  -  (ku  e15  sin  0   -  "Ju     cos   p) 


Ctf) 


From  (48)   and  (49)   we  find  that  the  time  averaged  electromagnetic   energy 
density  is  given  by 

J(c|E|2  +  ^H|2)~   |(|Uoo|2co82p  +    |uJ2sin2P)  (50) 

Upon  inserting  (26)    into  (50)  we  obtain 

J(€|E|2  +  n|H|2)~  f     W  -   (Wkx)"1/5[k(a-y)   +  (ka)1/5    (cr(oo)   +  C^oo))     cos20 


(ka)1/3    (CR(0)   +Cl(o])     sin2p]l  (5l) 


At  the  geometrical  shadow  boundary  y=a  the  right  side  of  (51)  approaches 
e/8  as  ka  becomes  infinite.   Let  us  define  the  shadow  boundary  for  finite  values 
of  ka  as  the  place  where  the  energy  density  has  this  same  value  (e/8).     Then 
(51)  shows  that  the  shadow  boundary  is  the  straight  line  y=a+s   .   Here  s  ,   the 
shift  in  the  electromagnetic  case,  is  given  by 

SEM  =  ^CR(a>)  +  cI(°°)]cos2P  +  [CR(°)  +  CjCoJJsin^l  (ka)"2/5a 


cos  0  +  s  sin  0  (52) 


(27).   Since  CT  =  J$   CL  for  Z=0  and  Z=oo,  (52)  can  be  written  &i 
1       R 


L8 


sm   =  [-1.180574  cos2p  +  1.560771  sin2p](ka)"2/5a  (53) 

Just  as  In  the  scalar  case,  we  conjecture  that  the  shift  along  each  ray  of 
the  boundary  of  the  shadow  cast  by  an  arbitrary  object  Is  given  by  (52). 
The  radius  a  is  defined  as  in  the  scalar  case.  The  angle  p  is  given  in  terms 
of  the  tangential  and  normal  components,  E    and  E   ,  of  the  incident  field 
at  the  point  of  tangency  of  the  ray  and  the  surface  of  the  object,  by 

tanp=  |E^nC|/|EjnC|  (5M 

This  conjecture  can  be  shown  to  be  consistent  for  cylinders  of  arbitrary  cross 
section, as  in  the  scalar  case. 

The  relations  (48)  and  (49)  also  hold  far  from  the  cylinder  in  the  forward 

direction,  i.e.  in  the  direction  of  the  positive  x-axis.  This  is  so  because 

ikx     -l/2  ikr 
in  that  direction  u  ~  e    +  Br  '  e    where  B  is  some  constant,  and  r  ~   x. 

Thus  Vu  ~  ikui  and  this,  when  used  in  (46)  and  (47),  yields  (48)  and  (49). 

Upon  applying  the  electromagnetic  cross  section  theorem  to  (48)  and  (49)  we 

obtain  for  the  total  scattering  cross  section  in  the  electromagnetic  case 

(fmi   =  tf  cos2p  +  a   sin2p  (55) 

EM    00         o 

Here  6    and  cr   are  the  cross  sections  for  the  scalar  hard  and  soft  cases, 
o      00 

respectively.  Upon  using  (38)  for  a    and  6       in  (55)  we  obtain  the  result 


6m~k& 


1  +[CR(oo)cos2p  +  CR(0)sin2p](ka)~2/5 


(56) 


-4a     Jl  +[0.498077  cos2p  -  0.432120  sin2p](ka)"2/5 
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By  comparing  (56)  with  (52)  or  (53)  we  find  that  C   and  s_,  are 


related  by 


tfEM  ~  "f  +  (1+  ^  )_1sEm]  =  ^  +  5Al4102  SEM  (57) 

As  in  the  scalar  case,  this  relation  suggests  that  for  an  arbitrary  cylinder 
the  total  electromagnetic  scattering  cross  section  per  unit  length  is 


^  +  (1+  /5)'1[(bem)1+  (sm)2 


(58) 


v  +  [0.498077  cos2p  -  0.432120  sin2pjk~2/5(a^5  +  a^/5)  I 


two  radii  of  curvature  of  the  cylinder's  cross  sectional  boundary  at  its 
points  of  tangency  with  the  shadow  boundary.   Similarly  for  a  three  dimensional 
object  the  total  cross  section  should  be 

°EM~2      [a+CI+v7?)'1    6    b^t 

(59) 

A  +  k"2/5  [0.498077  &   a1/5cos2BdT  -  0.432120  0  a^sin^dx]  I 


Let  us  apply  (59)  to  an  object  of  revolution  with  a  plane  wave  incident 
along  the  axis  of  revolution.   Let  p  be  the  maximum  radius  of  the  object.  Then 
(59)  yields 


0™  ~2*p2  +  2It[CR(oo)  +  CR(0)X2/5al/5p 


(60) 
~2«p2  +  2jt[0.065957]k"2/5a1/5p 
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Upon  comparing  (60)  with  (^-3)   we  find  that  for  axial  incidence  on  an  object 
of  revolution 


m 


(a    +  d    ) 

v  o    co  ' 


For  a   sphere   P*=a  and   (6c)   becomes 


°EM~2,ta    I1  +  0>65! 


■  f 


This  result  agrees  with  the  known  result  [4,5]  ^   this  case,  which  is  a 
verification  of  our  conjecture. 
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